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The study of the broad-band emission of GHz-Peaked-Spectrum (GPS) radio galaxies is a powerful tool to investigate the 
physical processes taking place in the central, kpc-sized region of their active hosts, where the jets propagate and the lobes 
expand, interacting with the surrounding interstellar medium (ISM). We recently developed a new dynamical-radiative 
model to describe the evolution of the GPS phenomenon (Stawarz et al. 2008): as the relativistic jets propagate through 
the ISM, gradually engulfing narrow-line emitting gas clouds along their way, the electron population of the expanding 
lobes evolves, emitting synchrotron light, as well as inverse-Compton radiation via up-scattering of the photon fields from 
the host galaxy and its active nucleus. The model, which successfully reproduces the key features of the GPS radio sources 
as a class, provides a description of the evolution of their spectral energy distribution (SED) with the lobes' expansion, 
predicting significant and complex X-ray to 7-ray emission. We apply here the model to the broad-band SED's of a sample 
of known, X-ray emitting GPS galaxies, and show that: (i) the free-free absorption mechanism enables us to reproduce 
the radio continuum at frequencies below the turnover; (ii) the lobes' non-thermal, inverse-Compton emission can account 
for the observed X-ray spectra, providing a viable alternative to the thermal, accretion-dominated scenario. We also show 
that, in our sample, the relationship between the X-ray and radio hydrogen column densitities, Nh and Nhi , is suggestive 
of a positive correlation, which, if confirmed, would support the scenario of high-energy emitting lobes. 
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1 Introduction 

It is currently accepted that the GPS galaxies sample the 
youngest fraction of the population of powerful radio ga- 
laxies. From sub-kpc scales, their jet-lobe's structures pro- 
pagate through the host-galaxy ISM, evolving into sub-ga- 
lactic (Compact Steep Spectrum, CSS) sources, which then 
expand to super-galactic scales (see O'Dea 1998 for a re- 
view). However, this scenario still has several open issues, 
like the absorption mechanism responsible for the characte- 
ristic turnover in the radio spectrum, the details of the dy- 
namical evolution and interaction with the ISM, the parame- 
ters of the central engine, and the origin of the high-energy 
emission. We recently proposed a model which addresses 
some of these issues through the analysis of the broad-band 
emission of GPS galaxies (Stawarz et al. 2008). 

Here we show that our model appears to be promising, 
enabling us to reproduce a number of observed properties 
of a sample of X-ray emitting GPS galaxies. 
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2 The model 

We recall below the key features of our dynamical-radiative 
model, referring the reader to Stawarz et al. (2008), and re- 
ferences therein, for a more comprehensive discussion. 

Our description of the dynamical evolution of GPS sour- 
ces mounts on the model proposed by Begelman & Cioffi 
(1989) to explain the expansion of classical double sources 
in an ambient medium with density profile p{r). The rele- 
vant equations can be derived by assuming that: (i) the jet 
momentum flux (proportional to the jet kinetic power Lj) is 
balanced by the ram pressure of the ambient medium spread 
over an area Ah, (ii) the lobes' sideways expansion velo- 
city equals the speed of the shock driven by the overpres- 
sured cocoon, with internal pressure p, in the surrounding 
medium, and (iii) the energy Ljt transported by the jet pair 
during the source lifetime is converted into the cocoon's 
internal pressure. For young GPS sources with age t, li- 
near size LS{t) ^ 1 kpc, and transverse size lc{t), expan- 
ding in the central core of the gaseous halo of the ellipti- 
cal host galaxy, we could constrain the model with a num- 
ber of reasonable assumptions: (i) a constant ambient den- 
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sity p = mpUo (with rUp the proton mass, and no ~ 0.1 
cm~'^ ), representative of the King-profile's core, (ii) a con- 
stant hot-spot advance velocity t;h — 0.1c, as suggested by 
many obervations of compact symmetric objects (but see 
Kawakatu, Nagai, & Kino 2009 for alternative scenarios), 
(iii) a scaling law lc{t) ^ t^/^, reproducing the initial, bal- 
listic phase of the jet propagation. All the lobes' physical 
quantities become thus functions of two parameters only: 
the jet kinetic power Lj and the source linear size LS. 

We then studied how the broad-band radiative output 
of GPS sources evolves as the source expands, for a given 
jet power Lj. The magnetic field in the expanding lobes 

scales as B = {^tttibpY^'^ ~ Ly^LS'^/"^ , with r/s = 
Ub/p < 3, and Ub the magnetic energy density. The elec- 
tron population Q(j) (with 7 the electron's Lorentz factor), 
injected from the terminal jet shocks to the expanding lobes, 
evolves under the joint action of adiabatic and radiative en- 
ergy losses, yielding a lobes' electron population iVg(7), 
which has a broken power-law form with critical energy 
7ci when Q{j) is a power law, and a more complex form 
when (5(7) is a broken power law with intrinsical break 
7int — 2 10"^. Assuming that the lobes' electrons, in rough 
equipartition with the magnetic field and the cold protons, 
provide the bulk of the lobes' pressure, the electron energy 
density is Ug = rjeP, with rje ^ 3. The lobes' electrons are 
source of synchrotron radiation, with luminosity I/gyn con- 
stant with time, and energy density Ugyn ~ LS~^/'^. Free- 
free absorption (FFA) of this radiation by neutral-hydrogen 
clouds of the narrow-Une region (NLR), engulfed by the ex- 
panding lobes and photoionized on their surface by the ra- 
diation from the active nucleus (as proposed by Begelman 
1999), is the process which we favour for the formation of 
the inverted spectra. While the synchrotron-self-absorption 
(SSA) process does not enable us to reproduce the observed 
turnover frequencies Vp, the spectra below the turnover, and 
the fp - LS anticorrelation (O'Dea & Baum 1997), FFA 
effects best fit the inverted spectra, and are a promising can- 
didate to account for the above anticorrelation. 

The lobes' particles also produce inverse-Compton (IC) 
radiation via up- scattering of both the synchrotron radia- 
tion (synchrotron-self-Compton mechanism; SSC) and the 
local, thermal photon fields generated by the accretion disc, 
the torus, and the stellar population of the host galaxy. The 
energy density U^ad of the thermal fields was evaluated by 
assuming that the nuclei of GPS sources share the proper- 
ties of quasars and Seyfert galaxies. The accretion disc, as- 
sumed to produce the bulk of its luminosity (10'*'''~'*^erg 
s~^ ; e.g., Koratkar & Blaes 1999) at UV frequencies, pro- 
vides C/uv ~ LS~^; the dusty torus, radiating the disc's 
UV photons at IR frequencies with efficiency e ^ 10 — 
100%, yields Uir ~ LS~^; finally, the host galaxy con- 
tributes near-IR to optical photons with f/opt independent 
of LS. The IC scattering of all the above radiation fields 
yields significant and complex high-energy emission, from 
X-ray to 7-ray energies. Whereas in GPS quasars the direct 
X-ray emission of the accretion disc's hot corona and of 



the beamed relativistic jets may overcome the X-ray output 
of most of these sources, in GPS galaxies those contribu- 
tions are expected to be obscured by the torus and Doppler- 
hidden, respectively, and the lobes are expected to be the 
dominant X-ray source. 

3 Observational supports 

Our model, and specifically the prediction of the X-ray- 
emitting lobes, may be supported by further observational 
evidence, which we did not discuss in Stawarz et al. (2008). 

The X-ray emission of GPS galaxies has been tradi- 
tionally interpreted as thermal radiation from the accretion 
disc, absorbed by a gas component associated with the AGN 
and characterized by an equivalent hydrogen column den- 
sity Nh (O'Dea et al. 2000; Guainazzi et al. 2004, 2006; 
Vink et al. 2006; Siemiginowska et al. 2008), rather than as 
non-thermal emission from the jets or the lobes. The for- 
mer scenario is mainly based on the apparent discrepan- 
cy between the equivalent total-hydrogen column density 
Nh derived from the X-ray spectral analysis and the neutral- 
hydrogen column density Nhi derived from the 21 -cm ra- 
dio measurements. Because Nh always appeared to exceed 
Nhi of 1-2 orders of magnitudes, it came natural to interpret 
the X-rays as produced in a source region which is more ob- 
scured than the region where the bulk of the radio emission 
comes from, and thus located in between the radio lobes; 
otherwise, an unreasonably high fraction of ionized hydro- 
gen (H II) would be necessary to account for the above dif- 
ference (e.g., Guainazzi et al. 2006; Vink et al. 2006). Such 
a scenario would also be consistent with the observed anti- 
correlation between Nhi and linear size found by Pihlstrom, 
Conway & Vermeulen (2003), being the fraction of ionized 
gas likely low in a young radio source with still expanding 
Stromgren sphere (Vink et al. 2006). 

The discrepancies between the Nh and Nhi values men- 
tioned above should actually be regarded with caution. The 
Nhi estimate is derived, from the measurements of the H I 
absorption fines, as a function of the ratio between the spin 
temperature Tg of the gas and its covering factor Cf , repre- 
senting the fraction of the source covered by the H I screen 
(e.g., Gupta et al. 2006). The common assumption Tg/cf = 
100 K refers to the case of complete coverage (cf = 1) of 
the emitting source by a standard cold (Tk — 100 K) ISM 
cloud in thermal equilibrium, and thus with spin tempera- 
ture equal to the kinetic temperature (Tg = Tk). However, 
this assumption returns a value of Nhi which represents a 
lower limit to the actual neutral hydrogen column density 
(Pihlstrom et al. 2003; Vermeulen et al. 2003; Gupta et al. 
2006). In fact, in the AGN environment, illumination by 
X-ray radiation might easily raise Tk up to 10^ — 10^ K 
(Conway & Blanco 1995; Maloney, HoUenbach, & Tielens 
1996), making Tg raise accordingly (Davies & Cummings 
1975; Liszt 2001); a source covering factor smaller than 
unity would also increase the Tg/cf ratio. Both the above 
effects might lead to Nhi values fully consistent with the 
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Nh estimates. Finally, temperatures as high as several 10'' 
K would likely imply the presence of a non-negligible frac- 
tion of HII (Maloney et al. 1996; Vink et al. 2006), also 
contributing to relax possible residual column-density dis- 
crepancies. The consistency of Nh and Nhi would make the 
scenario of non-thermal X-ray-emitting lobes a viable al- 
ternative to the accretion-disc dominated model. 

Evidence is mounting that, in GPS and CSS sources, 
the H I absorption lines are not generated by a screen cove- 
ring the source uniformly: instead, they seem to originate in 
clouds of neutral hydrogen connected with the radio struc- 
tures of their jets and/or lobes, and possibly interacting with 
them (Morganti et al. 2004; Labiano et al. 2006; Vermeulen 
et al. 2006). The association of the bulk of the H I absor- 
ption with the optical emission lines currently supports the 
identification of the absorbers with the atomic cores of the 
NLR clouds, although the presence of H I elsewhere is not 
ruled out (Labiano et al. 2006; Vermeulen et al. 2006). 

In our GPS model, the NLR clouds are gradually en- 
gulfed by the expanding lobes: besides being responsible 
for the FFA of the radio photons, they might play an impor- 
tant role in the absorption of the lobes' X-ray radiation. 

4 Comparison with broad-band data 
4.1 ModeUing the SED's 

We tested our dynamical-radiative model on the 11 GPS 
galaxies currently known as X-ray emitters.^ In Fig. 1, we 
show, as an example, the modelUng of the intrinsic^ broad- 
band SED of B0108H-388, a source with LS = 41 pc. The 
SED data were derived from the hterature, and properly de- 
absorbed. The modelling of the complete SED sample will 
be presented elsewhere. 

The radio data of lERS B0108H-388 were modelled as 
synchrotron radiation produced by a lobes' electron popu- 
lation iVe (7) derived from the evolution of an injected hot- 
spots' population (5(7) ~ 7"*, with s = s\ = 1.8 for 
7 < 7int, and s = 82 = 3.2 for 7 > 7int. FFA effects en- 
able us to best fit the spectral behaviour at frequencies be- 
low the ~ 6 GHz turnover. The thermal emissions from the 
torus (IR), the disc (UV), and the host galaxy (optical-NIR) 
were modelled as black-body spectra with the appropriate 
frequency peaks (tiR = 0.5 10^^ Hz, tt,v = 2.45 10^^ 
Hz, and fopt = 2.0 10^* Hz) and bolometric luminosities 
(LiR = 5.0 10'"' erg s^^ , Luv = 5.0 10''''^ erg s"! , and 
-^opt = 6.0 10^^ erg s~^ ). The comptonization of the syn- 
chrotron and thermal radiation fields yields the high-energy 
spectral components. For this source, the X-ray emission is 
dominated by the comptonization of the IR radiation. Our 

1 lERS B0026+346, lERS B0108+388*, lERS B0500+019*, lERS 
B07104439, PKS B0941-080, lERS B1031+567*, lERS B1345+125*, 
rVS B1358+624*, lERS B1404+286, lERS B2128+048, lERS 
B23524495*. Sources marked with an asterisk are those of the sub- 
sample discussed in Sec. 4.2, and included in Fig. 2. 

^ Throughout this paper, we use the cosmological paramters: Ha = 
0.7, Hm = 0.3, with i/o = 72 km s'^ Mpc"! 
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Fig. 1 Modelling of the intrinsic SED of GPS galaxy lERS 
B0108+388. Radio to X-ray data were derived from: NED; Dal- 
lacasa et al. (2000); Tinti et al. (2005); Stickel et al. (1996); 
Vink et al. (2006). The curves show the modelled spectral com- 
ponents: synchrotron emission, and corresponding SSC emission 
(solid lines); thermal star light (dotted line); comptonized ther- 
mal emission from the torus and the disc, respectively (dashed 
and dash-dot-dotted lines). The comptonized starlight's limiinosity 
does not appear in the plot because it is lower than 10''^ erg . 

model well reproduces the observed X-ray spectrum, and 
predicts significant 7-ray emission. 

4.2 Nh -Nhi connection 

Besides the modelling of the broad-band SED's, a way of 
discriminating among different scenarios, and unveil the ac- 
tual X-ray production site, is to compare the properties of 
the X-ray and radio absorbers, i.e. the Nh and Nhi column 
densities. Such a comparison can be performed either for 
individual sources, where an ad hoc increase of the Tg para- 
meter can remove possible Nh and Nhi discrepancies, or for 
a source sample, where the existence of a positive, signifi- 
cant Nh -Nhi correlation would suggest that the X-ray and 
radio absorbers coincide, thus supporting the co-spatiahty 
of the X-ray and radio source. 

We investigated the existence of a connection between 
Nh and Nhi in our sample. For a positive correlation, we 
searched the source subsample for which both Nh and Nhi 
estimates (either detections or upper limits) are available 
(see Fig. 2, and footnote 1). We obtained the following re- 
sults: (1) the subsample of 5 sources for which both Nh and 
Nm detections are available displays a strong (Pearson's 
r = 0.997) and highly significant (S = 2.3 • 10"*) Nh - 
Nhi positive correlation; (2) in the above-mentioned sub- 
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Fig. 2 Nh vs. Nhi for 6 GPS's of our sample^. Solid symbols: 
Nhi was computed with Ts = 100 K; arrows are upper hmits. 
Open symbols: as an example, the same sources with Ts = 6 10^ 
K are shown. Dash-dotted line: linear fit to the 5-source subsample 
of Nh /Nhi detections (with Ts = 100 K); dotted line: linear fit to 
the 6-source subsample including both detections and upper limits; 
Data are from: Guainazzi et al. (2006); Mirabel (1989); O'Dea et 
al. (2000); Pihlstrom et al. (2003); Siemiginowska et al. (2008); 
Vermeulen et al. (2003); Vink et al. (2006). 

sample, the strength and significance of the correlation sub- 
stantially decrease when using non-parametric methods 
(Kendall's and Spearman's correlation coefficients; e.g., 
Press et al. 1992); (3) the 6-source subsample including 
the Nhi upper limit also shows, according to survival anal- 
ysis techniques (ASURV, Rev. 1.2; e.g.. La Valley, Isobe & 
Feigelson 1992), a positive correlation, however with lower 
strength (p ~ 0.6), and a significance varying in the range 
S = 0.17 — 0.33, depending on both the Nh value (when 
more than one is available) and the statistical method cho- 
sen for the analysis. Although the data are suggestive of a 
positive correlation, further measurements would definitely 
help to improve the statistics. 

5 Conclusions and future prospects 

Our dynamical-radiative model can reproduce the observed 
broad-band SED's of X-ray emitting GPS galaxies. The 
shape of the radio spectra at frequencies below the turno- 
ver is best fitted by assuming FFA effects as the dominant 
absorption mechanism, whereas the X-ray spectra can be 
ascribed to IC scattering of the thermal radiation fields (ac- 
cretion disc, torus, and host galaxy) off the lobes' electron 
population. Further observational support to the X-ray-lo- 
bes' scenario comes from the radio and X-ray hydrogen col- 
umn densities of a sample of X-ray GPS's: the data suggest 
a positive correlation, which, if confirmed, would point to- 
wards the co-spatiality of the radio and X-ray emission sites. 
Additional measurements, necessary to improve the statis- 
tics, are already planned. 
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